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ABSTRACT 

We study the dynamics of the boundary in two dimensional dilaton gravity 
coupled to A^ massless scalars. We rederive the boundary conditions of [|l| 
and in a way which makes the requirement of reparametrization invariance 
and the role of conformal anomaly explicit. We then study the semi classi- 
cal behaviour of the boundary in the N = 24 theory in the presence of an 
incoming matter wave with a constant energy flux spread over a finite inter- 
val. There is a critical value of the matter energy density below which the 
boundary is stable and all the matter is reflected back. For energy densities 
greater than this critical value there is a similar behaviour for small values of 
the total energy thrown in. However, when the total energy exceeds another 
critical value the boundary exhibits a runaway behaviour and the spacetime 
develops singularities and horizons. 



1 Introduction 

In a recent paper [l|], we have studied a model of two dimensional gravity and 
massless matter field with a boundary in space time |0, ^. The reason for 
introducing such a boundary in the model is many-fold. In the two dimen- 
sional reduction of four dimensional spherically symmetric space time such 
a boundary may be thought as the origin of the radial coordinate. Moreover 
semiclassical studies of dilaton gravity models in two dimensions inevitably 
lead to a boundary beyond which the theory becomes strongly coupled. In a 
ficucial coordinate system {u, v) the boundary may be chosen to be the line 
u = V. We will refer to this coordinate system as the "fixed boundary gauge" . 
In these coordinates the conformal factor of the metric and the dilaton are 
independent fields. However, u and v are not asymptotically minkowskian. 
Alternatively we can work in Kruskal coordinates which are simply related 
to the asymptotically minkowskian coordinates. In these coordinates the 
dilaton is equal to the liouville mode and the boundary is dynamical. Its 
world line is determined by the matter energy momentum tensor, through a 
boundary equation derived in |I|, 0. In the fixed boundary gauge the same 
equation determines the relationship between the dilaton and the liouville 
mode. 

In the full quantum theory the boundary undergoes quantum fluctuations 
which are determined by the boundary equation, which has to be now re- 
garded as an operator equation. (In the fixed boundary gauge, these are the 
dilaton and liouville fields which fluctuate.) The dilaton-gravity sector of 
the theory has no physical propagating degrees of freedom, but the degree of 
freedom corresponding to the boundary survives. Thus we get an idea about 
the strength of quantum gravity in such simple model by studying boundary 
fluctuations, as done in [^. 

In some sense the boundary acts as a model for the high energy, strong 
coupling physics of the problem. Just as in the Callan-Rubakov effect we 
would like to model the high energy physics by suitable boundary conditions. 
The boundary conditions must respect some minimum requirements, like 
conservation of charges which should be conserved. In general such minimal 
requirements may not be restrictive enough to enable us to get insight into 
the real problem. However, it is also common that minimal requirements 
turn out to be rather restrictive and allow one to derive "model indpendent" 
consequences. The boundary conditions which arise in the Callan-Rubakov 



effect is such an example [Q. 

In the present problem, the boundary conditions must respect conser- 
vation of energy momentum and the general covariance of the underlying 
theory. In a gauge fixed theory the latter means that the reparametriza- 
tion constraints must be obeyed (in a conformal gauge this is the same as 
maintaining conformal invariance). 

The most general boundary conditions which respect these requirements 
were derived in 0, [l|. The conditions of general covariance and energy- 
momentum conservation turn out to be rather restrictive. Assuming that 
the matter is reflected off the boundary separately, the boundary conditions 
in the dilaton gravity sector is characterized by a single parameter (3, which 
is related to the value of the dilaton field at the boundary. In []I| a "boundary 
equation" was derived starting from the fixed boundary gauge. The behavior 
of the boundary and the resulting space-time structure is similar for all non- 
zero values of /3, while /3 = (which are the boundary conditions of [^ 0]) 
is a rather singular limit. In fact the boundary conditions arise dynamically 
when one adds a boundary cosmological constant term proportional to (3 to 
the two dimensional action 0, §], which immediately explains why /5 = is 
rather special. Our model is in fact similar to the boundary conformal field 
theories considered recently in |^. 

In |I[] the semiclassical boundary equation was solved in the N = 24 
theory for an incident shock wave. For (3^0 there is a critical value of the 
strength of the shock wave, i.e. the total energy. Below this critical value 
the boundary has a stable timelike evolution and no black hole is formed. 
When the energy exceeds the critical value the boundary runs away with 
ever increasing acceleration, approaching the speed of light. At the same 
time a black hole is formed which subsequently evaporates. The stability of 
the boundary for small disturbances have been recently analyzed in ^ for 
arbitrary A^ and an intriguing proposal based on the similarity of the model 
with open string theory has been put forward for the N = oo theory in |^. 

Note that a limit on the strength of the shock wave does not set any 
limit on the energy density for an aribitrary pulse. The /3 = theory, which 
is essentially the RST model, can be solved for arbitrary distribution of 
incoming matter energy and a black hole is formed when the energy density 
exceeds a critical value. However, in this model the boundary becomes space- 
like whenever this critical value is exceeded. It is of interest to know the 
behavior of the /? 7^ model for more general incoming matter distribution 



and ask whether there is some critical value for the energy density as well. 

In this letter we first rederive the boundary condition and the boundary 
equation for general A^ in the Kruskal gauge. This derivation makes the role 
of reparametrization invariance and the anomaly in the boundary conditions 
explicit. We then solve the semiclassical boundary equation for A^ = 24 for an 
incoming matter distribution which consists of a constant energy density in 
a finite interval of retarded time. We find that for the boundary to become 
unstable and a black hole to form, there is a critical value of the energy 
density as well as the total energy. 

2 The Boundary Equation 

The classical action of the model is given by 

S=^j dudv^g[e-'t>{R - A{V<j>f - 4A^ + \iy ff] (1) 

where gab{u,v), (/){u,v) and f'^{u,v) are the two dimensional metric, dilaton 
and the matter fields. R is the curvature scalar and A is the comological 
constant which we set to one by properly choosing the scale. 

We will use the conformal gauge where gab = e'^^Vab- Since the matter is 
conformally coupled its quantum effects are entirely contained in the Weyl 
anomaly which means that there is a liouville term in the action. There 
are additional liouville terms coming from the measures of integration of the 
conformal mode p, the dilaton and the ghosts. These measures are chosen so 
as to ensure that the Hawking radiation consists of only physical propagating 
particles, which are the matter fields /*. Following RST one can use the 
freedom to add a local counterterm to make the semiclassical model solvable. 
The resulting action may be written as 

S=- f dudv[2ndud,p- e^P-2{l- -)dupd,p+lduf d J' + (ghosts)] (2) 
n J K 2 

where we have defined the fields 

Q = -[e-^'l^ + Kcj)] p = p-(p (3) 



K 



and /« = H- 



It is also possible to write down models which are exactly confornially 
invariant f^, 10|. These models have an action which is identical to (Q). 



However the expressions for f2 and p in terms of (j) and p are different. The 
boundary conditions and the boundary equation which we will derive are 
valid for these theories as well. However the interpretation of the solution in 
terms of the space-time structure would be different. 

The equations of motion following from (^ are very simple 

dAp = dud.f = dA^ = -e^P (4) 

Following we introduce two chiral fields X^{u) and X~{v) and solve the 
equation for p 

^^P-^^ = d^X+{u)d,X-{v) (5) 

The equation for Vt has the general solution 

VL = -X+{u)X-{v)+g+{u)+g^{v) + K (6) 

where i^ is a constant and g^ are arbitrary chiral functions. We can trade 
these functions for two other chiral fields Y^{v) and Y~{u) which are defined 

as 

d^g+{u) = Y-{u)duX+{u), d,g-{v) = Y+{v)d,X~{v), (7) 

In the gauge fixed theory the requirement of general covariance of the 
underlying model is imposed by setting the total energy momentum tensor 
to be zero 

rpg I rpm _|_ rpgh _ Q T^S _|_ T"™ _|_ T^S^ = f) ('fl\ 

where j'a ^xm^rpgh g^^g^^d for the graviton-dilaton, matter and the ghost en- 
ergy momentum tensors respectively. In what follows, the ghosts completely 
decouple. The matter part T"^ is standard, while the graviton-dilaton part 
has the following expressions 

Tl = t,duY-duX++'^[{du\ogduX+y-2dl\ogduX+] 

Tl = Kd,Y^d,X- + '^[{d,\ogd,X-)^-2d^Jogd,X~] (9) 

where k = ^^i^ ■ ^o^^ that in (H) above, the matter part T"* contains a vac- 
uum energy term coming from normal ordering effects. Alternatively in the 



semiclassical theory the vacuum has to be defined in terms of modes which 
are positive frequency with respect to asymptotically minkowskian coordi- 
nates and there is a term which comes from the anomalous transformation of 
the energy momentum tensor to some other globally valid coordinate system. 
We now derive the general boundary condition. The minimal require- 
ments which the theory with these boundary conditions must satisfy are 

1. The dilaton field must be a constant along the boundary 

2. The matter must be perfectly refiected from the boundary 

3. The reparametrization constraints ® must be obeyed. 

The last condition guarantees that the boundary conditions are consistent 
with the general covariance of the model. 

Consider a timelike boundary in the theory specified by the equation 

V = fB{u) (10) 

or its inverse equation 

u = gBiv) (11) 

The condition (2) above and the constraint (H) imply that the total energy 
momentum tensor must be perfectly refiected. This means that the right 
moving energy momentum tensor must be set to be equal to the left moving 
energy momentum tensor after performing a coordinate transformation by 
the function (pH]). In doing so one must remember that energy momentum 
tensors transform anomalously with a Schwarzian derivative term. 
We will impose Dirichlet boundary conditions on the matter fields 

f{gB{v),v)=0 (12) 

We will also impose similar boundary conditions on the ghosts. This means 
that general covariance requires that the gravity energy momentum tensor 
must obey 

Tliv) = id.gBiv)f T^^[gB{v)] + k{gBiv),v} (13) 

where {., .} denotes a Schwarzian derivative 

{f,^} = y-^ijr (14) 
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In the bulk theory independent reparametrizations of u and v are still 
symmetries in the conformal gauge. This is no longer true if we choose the 
boundary to be a fixed line u = v as in the fixed boundary gauge. However, 
one of the conformal reparametrizations still remain a symmetry and we 
can use this to choose X~^{u) = u. This leaves X^{v) as the dynamical 
degreee of freedom, which is determined by the constraints. This is the 
gauge used in |@, Q. In the semiclassical theory, the constraints determine 
the function X~ {v) in terms of the energy momentum tensor of the matter. In 
the quantum theory the fluctuations of X~{v) are determined in terms of the 
matter fluctuations via the operator constraints. Since in this gauge QBiy) = 
1 the Schwarzian derivative term vanishes. However the nontriviality of the 
transformation properties of the energy momentum tenaor are contained in 
the expression for T-j^ in (H). 

Alternatively one may fix the Kruskal gauge X'^{u) = u and X^{v) = v. 
The boundary cannot be fixed to be some predetermined line any more. 
Rather the boundary curve is now dynamically determined, as we shall see 
soon. In the following we will use the Kruskal gauge and accordingly replace 
the indices u and v hj + and — respectively. 

The condition that the dilaton is constant along the boundary follows 
from the solution for Q in (H) and (0) 

(r-(x+) - X-) + (a+/^)(y+(x-) - x+)U-=/,(x+) = o (is) 

Let us introduce a new quantity 



hB{X+) = V5+/b(X+) (16) 

The condition ([T5|) may be satisfied by setting 

y+(A'-) = jB(A-) + F|/iB(9fl(X-)),X-| 

where F[hB{x), x] is some general functional of hsix) and a function of x. 

We now use these expressions for Y^ to evaluate rf_,_,Tf_ in (|13D and 
obtain a functional differential equation for F. In the Kruskal gauge one 
simply has T^+ = d^Y^ etc. However since the boundary curve is nontrivial, 
the term in (0) which involves k is nonvanishing. One gets the equation 

d+F + d+hs^ = (a+ log hB)F + -[{0+ log hnf - dl log hs] (18) 



This has the most general solution 



K, 



F[hB{X+),X+] = (3hBiX+) + -a+ log hs (19) 

where f3 is an arbitrary parameter. This completes the derivation of the most 
general boundary condition which satsifies the requirements stated above. 
The term involving k is clearly the result of the anomalous transformation 
law of the energy momentum tensor. 

Finally, using the above expression for F one may derive the equation of 
motion of the boundary. This is simply the statement that the total T++ 
vanishes and is given by 

Khl + K(3d+hB + kdl log hB + T^+ - — ^ = (20) 

The last term is the standard vacuum energy which appears in the Kruskal 
gauge. 

So far we have not specified the value of the dilaton field on the boundary. 
At this stage this may appear as a second parameter specifying the boundary. 
However, as we shall se below, this is actually determined by the condition 
that the vacuum solution of the boundary equation corresponds to flat space 
with a linear dilaton. 

The boundary equation ( pOD determines the shape of the boundary curve 
in terms of the matter energy momentum tensor. Once this is known, we 
can use the expressions for Y~ , Y^ to write the solution for Vt and hence 
derive the space-time structure. In the following we will concentrate on the 
behaviour of the boundary itself. 

It is useful to rewrite this equation in terms of the quantity 

qB{X+)=X+hB{X+) (21) 

In terms of qB the equation ( pUD becomes 

-(X+9+)2 log qB H^qB - l){X+d+)\og qB 

As noted in this represents the equation of motion for a particle driven 
with some external force and subjected to some damping. In vacuum Tj^l = 



and the external force is the gradient of a potential. Consequently one may 
expect that the ground state is characterized by a constant qb which is the 
minimum of this potential. However, for this to occur for real Qb one must 
have /3^ + 2 > 4-. For general N one cannot, therefore, choose /? = [Q]. 
Furthermore for the vacuum solution to be stable under small deviations, 
one must have k > and 3 > ^^^. 

In the rest of the paper we shall limit ourselves to A^ = 24 for which the 
vacuum is perturbatively stable. 

3 Semiclassical Solution with a finite dura- 
tion wave 

We now investigate the stability of the boundary for arbitrary incoming en- 
ergy and energy density for A^ = 24. Even in this case, we have not been 
able to solve the equation (P^ ) for arbitrary r++. However to understand the 
nature of threshold of instability and black hole formation we only need to 
have a situation where we can tune the energy density and the total energy 
independently. We will thus solve the above equation for an incoming mat- 
ter wave which has a constant energy density for a finite interval of retarded 
time. Thus we have a matter energy momentum tensor of the form 

T^i = ^(^[^(^^ - 1 + e) - 0iX^ - 1 - e)] (23) 

which corresponds to a constant energy density ^ in the interval 1 — e < 
X+ < 1 + e. The total energy is thus M = ae. 

In the following we will name the various regions of the X+ space as 
follows (i) Region I : X+ < 1 - e (ii) Region IIl-e<X+<l + e (iii) 
Region III X+ > I + e. 

In Regions I and II the solution is given by 

,.(.Y^)4l/3 + v/A|f;|;±^| (24) 



where d is an integration constant. In (^^ we have made the following 
definitions 



A = /32 + 2 S=^ (25) 



We want to have initial conditions such that the space time is in the 
vacuum, and the boundary has a vacuum hyperbohc solution discussed above, 
before any matter has come in. This means that in region I we have to choose 
d = 0. This determines the value of the dilaton field on the boundary to be 
the quantity [0] 

e-^%undary = ^[l + P{P + V^)] (26) 

Thus the value of the dilaton field is also determined in terms of the parameter 
f3 as mentioned above. 

The nature of the solution in region II is qualitatively different for small 
and large a regardless of the value of e. Let us define ao — ^ 
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3.1 a < ao 

For a < ao the solution in region II is given by 



where we have defined 

T = A-2a ^ = ^ (28) 

The constant of integration du is determined by requiring that qs is contin- 
uous across X"*" = 1 — e. This gives 

dn = (1 - e)" ^^^^ (29) 

From the solution ( P?] ) it is clear that qb and hence fs will become infinite and 
the boundary will run away if {X^)^ — dij = for some value of X+ > (1 — e). 
This requires du > [1 — e)^. However, using the definitions of A and T it is 
clear from (p9|) that du is always less than (1 — e)*^. Thus there cannot be 
any runaway of the boundary in the region II. 

The solution in region III is given by ( p^ with the integration constant 
denoted by dm. This is in turn determined by the solution in region II by 
matching across X^ = 1 + e, which yields the following equation for dm 



l + ey + dm VT (1 + e)" + dn 



[1 + eY - dm VA (1 + eY - djj 
10 



Since dn < (1 - e)'' < (1 + e)'' it is clear from (^ that dm < (1 + e)^. The 
solution (1^) then shows that qs cannot blow up anywhere in region III as 
well. 

We thus conclude that for a < y the boundary is stable and remains 
timelike throughout. 

3.2 a > aQ 

For a > ckq the solution in Region II is given by 

qB = ]^[{3-^ tan(/i log X+ + Cn)] (31) 

where we have defined 



S = 2«-A ^=^ (32) 

The intgration constant c// is once again determined by matching the solution 
with the vacuum solution at X^ = 1 — e. This gives 

tan ^log(l - e) + 7 

tan cii = — — ^ (33) 

7 tan /ilog(l — e) — 1 

where we have defined 



The solution in region III is still given by (^) which we rewrite as 



Cm is determined by matching across X+ = 1 + e which yields 

(1 + 7^)tan[^ log(f±f )] 
CjH = (1 + ey .n,,_r...„ .i+,,i ^ (36) 



Let us first look for runaways of the boundary in the region III. The 
solution shows that this can happen if the value of cm is such that the 
denominator in (|35| ) vanishes for a value of X~^ > (1 + e). This means that 

11 



one must have cm > (1 + ^Y- The solution for cm then shows that this 
would require ,for 7 > 1, 

- -^ < tan[/i log(l±^)] < -- (37) 

7^ — 1 1 — e 7 

whereas for 7 < 1 one must have either 



,1 + e,^ ^ 27 
1-72 



tan[/i log(^ :)] > ^^ (38) 



or 

tan[/i log(^)] < -- (39) 

i — e 7 

It may be easily checked that in the limiting case of a shock wave this 
reproduces the limit on the total energy derived in [^. This limit corresponds 
to a ^ 00 and e ^ with the total mass M = ae kept fixed. Then (35) 
reduces to the condition M > /3-\/A. 

For finite values of a and e one has a more complex behaviour, since the 
various integration constants are periodic functions of a and e. Let us study 
the behaviour of the boundary for some given value of a > ao and (3 and 
the extent of the pulse e increasing from zero. It follows from (0)-(0) that 
the conditions for runaway in region III are satisfied for e lying in the range 
el^) < e < el^) where 



2/i 



ijji 1 — 7"= 



e^^^ = tanhi — [n7r-tan-^(-)]} (40) 

where n = 1 , 2 , ■ ■ ■ and 

5n = (n - 1)9(1 - ^) + n6{j - 1) (41) 

In the above equation the branch of the inverse tangent function has been 
chosen to be the interval [— f , f ]. It is easy to check that e^^^ < e^f^ < e^li 

The behaviour of the boundary may be now summarized as follows. For 
e < ei , qb is finite and the boundary is stable. For e^ < e < e^ the 
boundary runs away in region III. For e[ < e < 62 the quantity qb diverges 
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in region II and there is no divergence in region III. When eg < ^ < e2 a 
new divergence of qb appears in region III, and the pattern continues. 

In the physical problem a divergence of qb results in a runaway of the 
boundary. Just as in the case of the shock wave analyzed in 0] a singularity 
develops which is asymptotic to the boundary where it diverges and there 
is a resulting event horizon. Thus for the physical situation at hand, the 
critical value of e for which the boundary starts running away and a black 
hole is formed is the lowest value of e for which any of the conditions (p?!)- 
( |38D are satisfied, i.e. e^ = el • When e exceeds a higher value e^ the 
runaway occurs inside the region of the incident pulse. The other runaways 
for higher values of e are irrelevant because they occur in a region which is 
already behind the singularity and hence not contained in the semiclassical 
space-time. Profiles of the boundary for fixed value of a and different e (and 
hence different mass M) are shown in the Fig.l - 3. 

Note that we have kept a, (3 fixed in the above discussion. It may be seen 
from (|40|) that Ccr decreases monotonically from 1 at a = ao to zero for infinite 
a, with ecr — >■ "^^ as a — > cxd. Thus a limit of e for some given a means that 
there is a lower limit on the total mass of matter Mcr{(y) = aCcr which has 
been thrown in. For a = ao one has Mcr = ao- As a increases, Mcr{ct) first 
rises to a maximum value and then decreases, eventually becoming asympotic 
to the constant shock wave value /3\/A for large a. 

We thus conclude that for the boundary to run away and a black hole to 
form there is a minimum value of the incoming energy density as well as a 
minimum value of the total energy. 

4 Outlook 

The fact that the runaway of the boundary occurs at the same threshold 
at which a black hole singularity is formed in the semiclassical model is 
significant. While we have demonstrated this in the N = 24 model, we 
expect that a similar behaviour would hold for N ^ 24 as well. However the 
differential euqtion in the latter case is higher order and the stability analysis 
would be significantly different. It would be interesting to get an idea of this 
by using numerical integration techniques. Some work in this direction has 
been done in [0. 

For the shock wave solution we found in |]T| that the runaway of the 
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boundary in the semiclassical theory also signals the threshold beyond which 
an aypnitotic observer measures large quantum fluctuations of the boundary 
near the horizon. It would be interesting to understand the nature and 
implications of the boundary flcutuations in a more general setting. We 
leave that for future investigation. 
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FIGURES 



Fig.l: Typical boundary curve for a < aQ for any e. Even for a > aQ and 
M < Mcr{a) the boundary curve is similar. No singularities are formed. 

Fig. 2: Boundary curve for a > ao and M > Mcr{oi). Boundary runs away 
after the matter infall ends. Cosequently space-time singularity is formed. 

Fig. 3: For yet higher value of M, i.e. for the corresponding e > e^^-*, the 
boundary runaway occurs in the region of incident pulse. 
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